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ABSTRACT 

Context. Submillimeter spectroscopic observations of comets provide an important tool for understanding their chemical composition 
and enable a taxonomic classification. 

Aims. We aim to determine the production rates of several parent- and product volatiles and the '^C/'^'C isotopic carbon ratio in the 
long-period comet C/2004 Q2 (Machholz), which is likely to originate from the Oort Cloud. 

Methods. The line emission from several molecules in the coma was measured with high signal-to-noise ratio in January 2005 at 
heliocentric distance of 1.2 AU by means of high-resolution spectroscopic observations using the Submillimeter Telescope (SMT) at 
the Arizona Radio Observatory (ARO). 

Results. We have obtained production rates of several volatiles (CH3OH, HCN, H'^CN, HNC, HjCO, CO, and CS) by comparing the 
observed and simulated line-integrated intensities. We calculated the synthetic profiles using a radiative transfer code that includes 
collisions between neutrals and electrons, and the effects of radiative pumping of the fundamental vibrational levels by solar infrared 
radiation. Furthermore, multiline observations of the CH3OH J = 7-6 series allow us to estimate the rotational temperature using the 
rotation diagram technique. We find that the CH3OH population distribution of the levels sampled by these lines can be described by a 
rotational temperature of 40 ± 3 K. Derived mixing ratios relative to hydrogen cyanide are CO/CH30H/H2CO/CS/HNC/H''CN/HCN 
= 30.9/24.6/4.8/0.57/0.031/0.013/1 assuming a pointing offset of 8" due to the uncertain ephemeris at the time of the observations 
and the telescope pointing error. 

Conclusions. The measured relative molecular abundances in C/2004 Q2 (Machholz) are between low- to typical values of those 
obtained in Oort Cloud comets, suggesting that it has visited the inner solar system previously and undergone thermal processing. The 
HNC/HCN abundance ratio of ~ 3. 1% is comparable to that found in other comets, accounting for the dependence on the heliocentric 
distance, and could possibly be explained by ion-molecule chemical processes in the low-temperature atmosphere. From a tentative 
H"CN detection, the measured value of 97 ± 30 for the H'^CN/H'^^CN isotopologue pair is consistent with a telluric value. The 
outgassing variability observed in the HCN production rates over a period of two hours is consistent with the rotation of the nucleus 
derived using different observational techniques. 

Key words. Comets: individual: C/2004 Q2 (Machholz) - molecular processes - radiative transfer - submillimeter: planetary systems 
- techniques: spectroscopic 



1. Introduction 

Comets spend most of their lifetime in the outer solar system 
and therefore have not undergone much thermal processing. Line 
emission from cometary atmospheres at submillimeter and radio 
wavelengths is a very useful tool for studying their physical and 
chemical conditions and relation with other bodies in the solar 
system ( [Biver et al?]|2002[ [Bockelee-Morvan et aL]|2004[ ). The 



* Based on observations carried out with the 10-m Submillimeter 
Telescope at the Arizona Radio Observatory, Steward Observatory, 
Mount Graham, Arizona, USA. 

"I" Current address: Department of Astrophysical Sciences, Princeton 
University, NJ 08544, USA 



coma structure and expansion velocity can be derived by fitting 
the observed line shapes using a molecular excitation code. In 
addition, mixing ratios of volatiles such as CH3OH, CO and CS 
can be compared with observed chemical abundances in proto- 
planetary disks to improve our understanding of planet forma- 
tion processes. 

The composition of comets has been investigated in the last 
two decades to develop a classification based on abundances of 
primary chemical species that displays a great compositional di- 



versity (A'Hearn et al. 



1995 Mumma et al. 2003 Bockelee 



Morvan et al.| 2004 1. More than 20 parent volatile species that 
release directly from ices in the nucleus, in addition to sev- 
eral radicals and ions formed by photodissociation in the coma. 
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have been detected via ground-based spectroscopic surveys at in- 
frared and submillimeter wavelengths and in situ measurements. 
The composition of some cometary ices show strong evidence 
of processing in the solar nebula and can provide clues about 
their place of formation and subsequent evolution. The abun- 
dance of HCN relative to water has been observed to be roughly 
constant with a value of 0. 1 % in several comets for a wide range 
of heliocentric distances (B iver et al.|2002 ). Other species show 
a wide spread of production rates. For instance, CH3OH has 
been found to have a variable abundance relative to water, rang- 
ing between less than 0.15% to 6% at different heliocentric dis- 
tances (Bock elee-Morvan et al.| [2004). There is no evident cor- 
relation between the observed relative abundances and the dy- 
namical class of the comets. Water is the most dominant volatile 
species and is typically used to determine relative abundances. 
Although water is not directly accessible from the ground at sub- 
millimeter wavelengths, it has been observed from space using 
the Submillimeter Wave Astronomical Satellite (SWAS), Odin 
and Herschel satellites (Neufeld et aL]|2000t [Lecacheux et al.| 
2003 |Hartogh et al. 2009") or inferred by observations of the hy 



droxyl (OH) radical at radio wavelengths (jCrovisieretaLp002|. 
HerscheVs, Heterodyne Instrument for the Far Infrared (HliFl) 
is able to determine water production rates accurately (Hartogh, 
et al. 2010| l. Direct measurements of water can also be performed 



using ground-based telescopes by observations of non-resonance 



fluorescence emission at infrared wavelengths (Mumma et al. 
[T986l|Bonevetal.|2009D . 

Comet C/2004 Q2 (Machholz) was discovered on 27 August 
2004 by Donald E. Machholz ( [Machholz et al.|2 004). The comet 
passed perihelion on 24 January 2005 at a heliocentric distance 
of rh = 1.205 AU and geocentric distance of A = 0.435 AU. 
During its closest approach to Earth at a distance of 0.35 AU on 
6 January 2005 it reached a naked-eye visual magnitude of ~ 
3.5 as reported in the International Comet Quarterly. Owing 
to its favorable viewing geometry in the northern hemisphere, 
the comet was extensively observed from the ground at various 
wavelengths. C/2004 Q2 (Machholz) is a long-period highly ec- 
centric comet whose origin is most likely the Oort Cloud ac- 
cording to the comet classification scheme by |Levison| ([T996|, 
with an approximate initial orbital period of 118 000 years and 
eccentricity of 0.9994658 before the comet was perturbed grav- 



itationally in the inner solar system (Marsden 2004 Nakano 



2006| l. These values are close to the perihelion osculating ele- 
ments (see jM arsden 2005, JPL Small-Body Database^. Its dy- 
namical classification is still a matter of debate because of the 
strong non-gravitational forces that make a backward orbital in- 
tegration very unreliable. 

In this paper we present high-resolution spectroscopic obser- 
vations of several volatiles from comet C/2004 Q2 (Machholz) 
acquired at the Submillimeter Telescope (SMT). Seven species 
are detected, namely CH3OH, HCN, HNC, 



H2CO, 



CO, CS, 



and a marginal detection of H'^CN. The comet was observed 
shortly pre-perihelion in January 2005 when it was at a distance 
~ 0.36 AU from Earth. These observations provide information 
about the outgassing of several molecules and an isotopologue 
of HCN relative to HCN, which is often used as a proxy for 
water in cometary taxonomies. We calculate the CH3OH rota- 
tional temperature of the ground vibrational level from several 
rotational lines and production rates for the observed molecules 
using a radiative transfer code to fit the observed line intensities. 
Section |2] presents our SMT observations of comet C/2004 Q2 
(Machholz) and the reduction method. In Sect. [3] the radiative 



transfer models and analysis of the observations are described. 
Finally, we discuss the obtained results in Sect. |4] 



2. Observations 

A spectral line survey of primary volatile species in comet 
C/2004 Q2 (Machholz) was made using the SMT telescope lo- 
cated at the Mount Graham International Observatory (MGIO), 
a division of Steward Observatory on Mount Graham, Arizona 
dBaars & Martini [T9961 |Baars et al.||T999] l. The SMT has a 
parabolic 10-m primary dish and a hyperbolic secondary reflec- 
tor The observations were performed with the 0.8 mm double- 



sideband receiver using various spectrometers ( Villanueva et al. 
2005 de Val-Borro et al. 201 1 1. Detected emission lines may 
have a different gain response depending upon which side- 
band the lines were observed in. Thus, the sideband gain ra- 
tio (Gusb/Glsb) deviates from unity and introduces an addi- 
tional uncertainty of ~ 10% in the absolute brightness temper- 
ature calibration. A typical system temperature of 150 K was 
attained during the observations. We observed the comet simul- 
taneously with the chirp transform spectrometer (CTS) with a 
bandwidth of 215 MHz, and the acousto-optical spectrometers 
(AOSA, AOSB and AOSC) wifli total bandwidths of 1 GHz, 
970 MHz and 250 MHz, respectively. The spectral resolution 
provided was 40 kHz for the CTS, and 934, 913 and 250 kHz for 
the AOSA, AOSB and AOSC. One of the purposes of the obser- 
vations was to test the performance of the newly installed high- 
resolution CTS built at the Max Planck Institute for Solar System 
Research ( Hartogh & Hartmann |1990[ Villanueva & Hartogh 
2004; Vi llanueva et al.|2006 1. High spectral resolution is crucial 
for resolving the shape of rotational lines in comets and study 
the gas velocity and asymmetries related to non-isotropic out- 
gassing. 

Comet C/2004 Q2 (Machholz) was observed near perihelion 
during six nights in the period 13-18 UT January 2005. Here we 
focus on the observations acquired on the first four nights with 
good observing conditions. The data were taken using the stan- 
dard position switching observing mode, where a reference sky 
position separated from the comet by 0°.5 was observed for the 
same amount of time and subtracted from the on-source observa- 
tions. We used a 30-second exposure for the source and reference 
positions with 8-minute scans. A summary of the observing log 
with the total integration times for each line is shown in Table[T] 

The main beam equivalent brightness temperature scale, 
r^B, was corrected for the beam efficiency of the telescope es- 
timated from observations of Mars and Saturn in the interval 
fall 2007 -spring 200^ a nd calibrated using the chopper-wheel 
method (Ulich & Haas 19761. The frequency scale was con- 
verted into Doppler velocities in the nucleus frame and corrected 
for the relative motion of the comet with respect to the observer. 
We used orbital elements provided by the JPL HORIZONS sys- 
terrj^ to track the comet and the ephemeris to calculate the po- 
sition and relative motion of the comet with respect to the tele- 
scope. Our observing strategy included pointing and calibration 
observations approximately every hour The pointing of the tele- 
scope was checked and corrected by dedicated reference obser- 
vations of bright sources like Saturn, the protobinary system 
W3(OH) and the massive star-forming region DR21(0H) be- 
cause they were close to the comet. Typical pointing errors were 



,http://kpl2m . as . ar izona . edu/ smt_docs/smt_beam_eff . 



http : //ssd. jpl .nasa. gov/sbdb . cgi?sstr=C/2884+Q2 



htm 



^ http : //ssd. jpl -nasa . gov/?horizons 
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Table 1, Log of the SMT obsei-vations of comet C/2004 Q2 (Machholz) in January 2005. 



Date" 


Molecule 


Transition 




Integration'" 


('•h>" 


(A)" 




Beam size 


(UT) 








(min) 


(AU) 


(AU) 






12.96 


f HCN 
H2CO 
CO 


4-3 
4-3 
3-2 


13 


107.2 


1.220 


0.361 


42?5 


2a'93 
21'.' 09 
21'.'46 


13.20 ■ 


CS 
H'^CN 
CH,OH 


7-6 
4-3 
13i-13oA- 


33 


270.6 


1.219 


0.362 


42!6 


21'.' 64 
21'.' 49 
21'.' 65 


13.95 


CH,OH 


7-6 


16 


133.0 


1.217 


0.365 


43° 3 


21'.' 93 


14.17 


HNC 


4-3 


22 


183.0 


1.217 


0.366 


43° 5 


2a'46 


15.14 


HNC 
CO 


4-3 
3-2 


39 


329.2 


1.215 


0.371 


44° 2 


2a'46 
21'.'46 


16.16 


CS 
H'^CN 
CH,OH 


7-6 
4-3 
13i-13oA- 


47 


378.0 


1.213 


0.376 


45° 


21'.' 64 
21'.' 49 
21'.' 65 



Notes. In some receiver tunings several lines were observed simultaneously in the lower and upper sidebands. Start times in universal time (UT) 
fractional day values. Number of individual scans. Total integration time. Heliocentric distance. Geocentric distance. Solar phase 
angle (Sun-C/2004 Q2 (Machholz)-Earth). 



~ 4", i.e., about a quarter of the half-power beam width (HPBW) 
at the observed frequencies. 

These observations were very challenging owing to the rela- 
tively uncertain ephemeris of C/2004 Q2 (Machholz) since they 
were carried out about five months after the comet was discov- 
ered, and particularly because of the strong non-gravitational 
forces needed to reproduce the observations afterwards. Its or- 
bit changed significantly during its passage through the inner 
solar system due to planetary perturbations from an initial or- 
bital parameter of a^i' - 0.00 0404 AU ' to a future value of 
0.001856 AU"^ ( |Marsden|2004) l. At the time of the observations, 
the comet was moving with a fast apparent motion of ~ 300" 
per hour in the sky relative to the background stars. This intro- 
duced an additional term of uncertainty in the pointing accuracy. 
Comparing the ephemeris actually used for the tracking with the 
HORIZONS ephemeris calculated from the after-the-fact orbit 
including non-gravitational accelerations, we obtain an average 
pointing offset of 4". Including the effect of the telescope point- 
ing error, the offset was at maximum 8". This value was consid- 
ered in our computation of the production rates for the observed 
species shown in Table |2] 

The data analysis was performed using the CLASS analy- 
sis software, which is part of GILDAS package]^ in conjunction 
with the NumPy and SciPy libraries of high-level mathematical 
tools (Jones et al. 2001 OIiphant"2007V A standing wave ap- 
pears as a baseline ripple in some of the spectra. We determined 
the baseline by fitting a polynomial to the emission-free back- 
ground and subtracted it from the original spectrum. After re- 
moval of the standing wave, the individual scans were averaged 
to increase the signal-to-noise ratio (S/N). Line intensities were 
calculated from the weighted averages of the spectra where the 
statistical weights are the inverse square of the root mean square 
(rms) noise of each individual spectrum. 

The comet's heliocentric distance was in the range - 
1.213-1.220 AU, and the geocentric distance was A = 0.361- 
0.376 AU during the observing period (see Table[T]l. The HPBW 
of the telescope varied between 20'.'5 and 22" at different fre- 
quencies, corresponding to 6 600 and 6 900 km projected on the 
comet, so that molecules from the outer coma contribute to the 
detected emission. The line excitation in this region is dominated 
by colUsions between neutrals and electrons, and infrared fluo- 



rescence by solar radiation. During the observing run the solar 
phase angle (Sun-C/2004 Q2 (Machholz)-Earth) ranged from 
42°5 to 45°0. We detected seven species (CH3OH, HCN, H'^CN, 
HNC, H,CO, CO, and CS), whose spectra are shown in Figs.fT]- 
[2] and [4]^ The observations were focused on the study of tne 
hydrogen cyanide chemistry with the detection of HCN, HNC 
and a tentative detection of H'^CN to determine the isotopic 
'^C/'"'C ratio, as well as other volatile species that sublimate di- 
rectly from the nucleus into the coma and abundant daughter 
species. 

In Table [2] we show a comparison of the mixing ratios 
of the detected species with respect to HjO with those de- 
rived from infrared observations on similar dates (Bonev et al.l 



http : //www . iram . f r/IRAMFR/GILDAS 



2009} IKobayashi & Kawakita||2009| . For species that are de- 
tected in several nights during our observing run, the mixing 
ratios refer to the weighted average of the observed spectra. 
Table |3] shows the molecules and transitions detected in our 
survey with 1-cr uncertainties. Line frequencies were obtained 
from the latest onUne edition of the JPL Molecular Spectroscopy 
Catalog ( [Pickett et al. 1998). Most transitions were observed 
with the higher resolution CTS, which resolves the line shape, 
and the AOS backends simultaneously. Integrated line intensi- 
ties shown in Table [3] were obtained integrating over velocities 
in the interval [-2, 2] km s 'unless another line overlaps within 
that region. Doppler shifts in the comet rest frame were ob- 
tained as the first moment of the velocity over the same inter- 
val (2, rmB,v',7 Yji TmBi whcrc V is the Doppler velocity and the 
index / refers to the channel number). 

3. Results 

3.1. Radiative transfer modeling 

We adopted a molecular excitation and radiative transfer model 
based on the Sobolev escape probability method, which includes 
collisional effects and infrared fluorescence by solar radiation to 
derive the production rates (Bockelee-Morvan 1987 ; Biver 1997J 
Biver et al.|1999 1. In the outer coma, solar infrared pumping of 
vibrational bands followed by instantaneous spontaneous decay 
establishes a fluorescence equilibrium. To test the effect of the 
infrared excitation, an independent molecular excitation model 
based on the publicly available accelerated Monte Carlo radia- 
tive transfer code ratran was used to calculate the populations 
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Table 2. Compaiison of production rates of detected species relative to HCN and HjO in comet C/2004 Q2 (Machholz) with 
statistical uncertainties. 



Molecule 




Q/QncN 




e/gH^o^ 






(molec. s~^) 




(%) 


(%) 


(%) 


HCN 


(Z.Zo ± U.U2) X U) 


1 


r\ r\o A 1 A A A 1 
(J.(Jo4 ± U.UUl 


" -^-0.02 


0.16 ± U.Ul 




d n -\- n 1 






J.U / ± U.J 1 




CH3OH 


(5.5 + 0.6) X 10" 


24.6 ± 2.5 


2.1 ±0.2 1 


2.14 ±0.12 
1.54 ±0.07^ 


1.2 ±0.1 
1.65 ±0.09' 


H2C0 


f (3.87 ± 0.20) X 10^* 
(1.09 ± 0.06) X 10" ' 


1.7 ±0.1 

4.8 ±0.3 


0.14 ±0.01 
0.41 ± 0.02 


0.11 ±0.03 


0.18 ±0.01 


cs 


(1.15 ± 0.04) X 10^^ 
(1.29 + 0.04) X 10^*^ 


(5.1 ± 1.7) X 10-' 
(5.7 ± 1.9) X 10-' 


(4.3 ± 0.1) X 10-2 
(4.8 ± 0.2) X 10-2 






HNC 


(7 ± 2) X lO^" 


(3.1 ±0.9) X 10-- 


(2.6 ± 0.7) X 10-^^ 






h"cn 


(3 + 1) X lO^" 


(1.3 ± 0.5) X 10-- 


(1.1 ±0.4) X 10-^^ 







Notes. Weighted mean of production rates measured at the SMT on 13-16 January 2005 derived from excitation and radiative transfer models 
assuming a pointing offset of 8". A water production rate of Qh^o = 2.7 x 10^' molec. s-' was used to derive the mixing ratios, intermediate 
between the measurement on 20 January 2005 by the Odin satellite jBiver et al.||2007b^ and the value derived from observations at infrared 
wavelengths on 19 January 2005 (Bonev et al. 2p09|. Weighted m ean of mixing ratios obtained on 28-29 November 2004 and 19 January 2005 
with NIRSPEC on the Keck II 10-m telescope |Bonev et al.|2009| ''^ Mixing ratios o btained on 30 January 2 005 with NIRSPEC on the Keck II 
10-m telescope 1 Kobayashi & Kawakita 2009 1 Revised CH3OH mixing ratios from 
CS assume the presence of an extended source described in the text. 



Villanueva et al. 



1 2012 1. Production rates for H,CO and 



of the rotational levels as a function of the distance from the nu- 
cleus for various molecules and line emi ssion in the cometary 
coma IHo gerheijde & van der Tak|[2000]p ] This code includes 



collisional effects with water molecules and electrons, but ne- 
glects the pumping by solar infrared radiation from the ground- 
state vibrational level. We used the one-dimensional spherically 
symmetric version of the code following the description outlined 
in|Bensch & Bergin] (|2004[|, which has been extensively tested 



and used to interpret Herschel cometary observations (see e.g. 
Hartogh et aLl[2UTol [de Val-Borro e'rar][20Tol [Hartogh et al 
2011| l7 

The radial gas density profiles for parent molecules HjO, 
CH3OH, HCN, H'^CN, HNC, and CO were obtained using 
the standard Haser spherically symmetric distribution for parent 
volatiles ( |Haser|1957] l: 



np(r) = exp 

47rr^Vexp \ Vexp 



(1) 



where Q is the total production rate in molecules s-'. The pho- 
todissociation rate /3p takes into account the dissociation and ion- 
ization of molecules by the solar UV radiation, Vexp is the ex- 
pansion velocity and r is the nucleocentric distance. However, 
second-generation molecules such as HNC, CS and HjCO are 
not correctly described by this profile. In the last two cases, the 
density profile of daughter species originating from the photodis- 
sociation of an extended parent source in the coma is given by 
( jCombi et al.[2004l i 



Q 

47rr2v(.xp j6d - y8p 



exp exp 



(2) 



where y6p and Pa are the parent and daughter photodissociation 
rates. We obtained the photodissociative lifetimes for CH3OH, 
HCN, HNC, HjCO and CO from |Crovisier| ( [T994| assuming the 
quiet-Sun reference spectrum, arid the CS radical lifetime was 
taken from Biver et al. (201 1). The photodissociation rates were 



scaled by the heUocentric distance of C/2004 Q2 (Machholz) 



^ the source code is available from ratmns website at |http : //mm . 
sron.rug.nl/-vdtak/ratran/frames.html| 



at the time of the observations. Deviations from a spherically 
symmetric Haser distribution due to anisotropic outgassing were 
not considered. 

Neutrals in cometary atmospheres are excited by collisions 
with other molecules and electrons, which are the dominant ef- 
fects in the inner coma, and radiative pumping of the fundamen- 
tal vibrational levels by the solar infrared flux. Infrared pumping 
of vibrational bands by solar radiation contributes to the excita- 
tion in the outer coma where the gas and electron densities are 
low (Bockelee-Morvan 1987). Most of the detected emission in 
comet C/2004 Q2 (Machholz) originates from a region in which 
molecules are in an excitation state intermediate between fluo- 
rescence and collision-dominated equilibrium. 

Collision rates with water and electrons were obtained for 
several molecules from the cuiTent v ersion of the Leiden Atomic 
and Molecular Databasf](LAMDA; |Sch6ier et al.|2d05] l. When 
collision rates with water were not available, we scaled the 
collision rates with molecular hydrogen by the ratio of their 
molecular weights. However, the calculated line intensities of 
the rotational transitions in the ground-based vibrational level 
are weakly dependent on the collision cross-sections - of about 
a few percent. 

We used a water production rate of 2h,o = 2.7 x 
10^^ molec. s-', close to the Odin satellite measurement on 20 
January 2005 of (2.64 + 0.08) x lO^** molec. s-' derived from 
the observation of the fundamental transition at 557 GHz ( Biver] 
|et al.|[2007bi l at - 1.208 AU, and to the water production 
rates of (2.73 + 0.07) x 10^^ molec. s-' and (2.76 ± 0.08) x 
10^^ molec. S-' derived from observations at infrared wave- 
lengths on 19 January 2005 ( |Bonev et al.][2006l [2007] [2009] l. 
Excitation parameters in the model are the neutral gas kinetic 
temperature, which controls the molecular excitation in the col- 
lisional region, and the electron density. We assumed a gas ki- 
netic temperature of 60 K and the rotational temperature de- 
rived from the relative line intensities of the / = 7-6 transitions 
for C H,OH. The electron density and temperature profiles from 



Biver (1997i were assumed. Since the electron density in the 



coma is not a fully known quantity, an electron density scaling 



http : / /www . strw . leidenuni v . nl/ -moldata/ 
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factor of = 0.5 with respect to the reference profile deduced 
from in situ measurements of comet IP/Halley was used. The 
expansion velocity is assumed to be constant in the coma, with a 
value of 0.75 kms obtained from the width of the HCN line. 
This value is consistent with the shape of several lines observed 
with the Institut de Radioastronomie Millimetrique (IRAM) 30- 
m radio telescope at about the same time (Biver et al. in prepa- 
ration). 

The radiative transfer equation was solved by integrating 
along various lines of sight through the coma that covered 
2.5xHPBW at a given offset from the nucleus, to obtain the 
brightness distribution in the plane of the sky once the level pop- 
ulations were calculated. The beam-averaged emission was then 
calculated at the distance of the comet. 

Production rates were calculated comparing observed line- 
integrated intensities with the integrated synthetic lines obtained 
with the Sobolev escape probability and ratran codes by fitting 
the production rate to the observed line intensity to derive the 
relative molecular abundances in the nucleus ( Bensch & Bergin 
|2004 l). Opacity effects are negligible for all lines, therefore the 
difference in the derived production rates in our models arise 
only from including the pumping of excited vibrational bands by 
infrared solar radiation as an excitation mechanism ( |Bockelee^ 
Morvan & Crovisier||1989 1. Table [2] summarizes the mixing ra- 
tios relative to water and hydrogen cyanide in comet C/2004 Q2 
(Machholz). We show the derived production rates with statis- 
tical uncertainties for the observed molecules using parent- and 
daughter density distributions in Table[3] Production rates are af- 
fected by pointing errors of up to 8" since the ephemeris used at 
the time of the observations was relatively uncertain. Therefore 
we considered this pointing offset to compute the production 
rates. 



3.2. CH3OH 

Methanol has been observed in several comets with a wide range 
of abundances since its initial discovery in comet C/1989 XI 
(Austin) by Bockelee-Morvan et al. ( |1990 1. Methanol rotational 
lines often appear in multiplets at millimeter and sub-millimeter 
wavelengths, allowing the estimation of the temperature and ex- 
citation conditions in the coma. Rotational states are ordered in 
non-degenerate A levels and degenerate E levels corresponding 
to different symmetry states. Radiative transitions between A 
and E species are strictly forbidden, and collision-induced tran- 
sitions are highly improbable. We detected twelve A^, A' and 
E-methanol emission lines of the J - 7-6 series in C/2004 Q2 
(Machholz) at frequencies listed in Table [3] In addition, there is 
a blend of three emission lines at 338.513 GHz and a blend of 
two lines at 338.722 GHz. Figure [T] shows the observed CH3OH 
rotational spectrum between 338. ID and 338.75 GHz with labels 
indicating the transitions listed in Table [3] Several unidentified 
emission features are found at about 338.178, 338.232, 338.277, 
and 338.651 GHz in the lower sideband and 341.169, 341.543, 
341.588, and 341.642 GHz in the upper sideband, which are la- 
beled by (?) in Fig[T] The line widths of these features are com- 
parable with the detected CH3OH lines and thus are consistent 
with a cometary origin, although no plausible lines were found 
at those frequencies nor the corresponding image band frequen- 
cies in the JPL Molecular Spectroscopy Catalog (Pickett et al. 
1998[l and the Cologne Database for Molecular Spectroscopy 
( [Muller et"aLl|200H |2005| l. Integrated intensities and velocity 



emission in the direction toward the observer, which should be 
studied with a non-spherically symmetric model. Additionally, 
we detected serendipitiously a high-excitation-level methanol 
transition J = 13i-13oA at 342.730 GHz in the CS/CO spec- 
tra from 13 and 16 January. Figure |2] shows the averaged CTS 
spectrum from both dates. 



3.2.1. Rotational temperature 

The rotational temperature was calculated from the relative in- 
tensities of the individual CH3OH lines between levels with 
quantum numbers J - 7-6 and 13i-13oA^ using the rotation dia- 
gram technique, assuming that the population distribution of the 
levels sampled by the emission lines is in local thermodynamical 
equilibrium (LTE), i.e., described by a Maxwell-Boltzmann dis- 
tribution characterized by a single temperature. Then the column 
density of the upper level within the beam, A^u, can be expressed 
as 

A^,=A^-|:^exp(-^i (3) 



z(r™t) 



7^1-0 



where is the degeneracy of the upper level, Z denotes the par- 
tition function, which is a function of temperature, Tiot is the 
rotational temperature, is the energy of upper state, repre- 
sents the Boltzmann constant, and is the total column density 
averaged over the beam. This method is commonly used in stud- 
ies of the interstellar medium. 

A linear fit of \n{Nu/gu) versus provides the rotational 
temperature as the inverse of the slope and the averaged total 
column density divided by the partition function from its inter- 
cept at Eti = 0. We included in the analysis the blended transi- 
tions at 338.722 and 338.723 GHz starting from the 72 and 7_2 
E levels, and the 338.541 and 338.543 GHz pair from the 73 A+ 
and 73 A" levels, by assigning half of the observed intensity to a 
virtual level of intermediate energy (see Bockelee-Morvan et al. 
1994bf The high-excitation-level line at 342.730 GHz was in- 
cluded in the rotational diagram analysis. For optically thin con- 
ditions, which normally apply to cometary lines of most volatile 
species, the observed intensity, J T^^ dv, is proportional to the 
column density of the upper transition level ( [Bockelee-Morvan 
|etal.|1994a| i: 



TmB dv. 



(4) 



shifts of the identified CH3OH lines with statistical uncertain- 
ties are shown in Table [3] The observed velocity shifts toward 
the blue wing suggest asymmetric outgassing with preferential 



where Aui is the Einstein coefficient for spontaneous emission 
and the transition frequency. 

The rotation diagram for multiple CH3OH lines that fits 
both A and E-CH3OH simultaneously was determined using 
a weighted linear least-squares method where the weights are 
equal to the reciprocal of the variance of each line intensity 
measurement. The standard deviation of the fitted parameters 
is obtained from the diagonal elements of the covariance ma- 
trix. Figure |3] shows the best fit for all observed lines and only 
the J - 7-6 series with statistical uncertainties. In contrast to 
water, no reliable A/E ratio was retrieved from radio lines for 
CH3OH ( Pardanaud et al.|2007 1. This ratio is almost insensitive 
to temperature, i.e., it sigmhcantly departs from the statistical 
ratio only for very low spin temperatures, because of the small 
energy difference of the A and E levels. Retrievals of the A/E ra- 
tio in CH3OH for comet C/2001 A2 LIN EAR at infrared wave - 
lengths are consistent with Tspin > 18 K ( Villanueva et al. 2012| l. 

The measured rotational temperature m C/20(j4 Q2 
(Machholz) is 85 ± 7 K with 1-cr uncertainty from the linear 
best fit of the rotational diagram that includes all the observed 
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Fig. 1. CH3OH averaged spectrum obtained on 13.95 UT January 2005. The lower and upper x-axis scales represent the frequency 
of the lower and upper sidebands, respectively. Labels indicate detected CH3OH spectral lines listed in Table [3] There is a blend of 
three and A emission lines at 338.513 GHz (label 5). The observed line at 338.722 GHz is a blend of two t components (label 
13). The A Hne at 341.416 GHz (label 14) is in the upper sideband close to the A+ transition at 338.409 GHz (label 3). A marginal 
detection of the E line at 338.530 GHz is indicated by label 6. The unidentified emission features labeled by (?) are discussed in 
Sect. 112] 




Fig. 2. Averaged spectrum of the CH3OH (13i-13oA ) line at 
342.730 GHz observed on 13.20 and 16.16 UT January with the 
CTS. The vertical axis is the calibrated main beam brightness 
temperature and the horizontal axis is the Doppler velocity in 
the comet rest frame. 



methanol lines on 13.20, 13.95 and 16. 16 UT January (solid line 
in Fig. [3]). The high-excitation line J = 13i-13oA" is too bright, 
indicating that the LTE assumption is not valid for that tran- 



sition. Weighted least-squares fitting methods are known to be 
very sensitive to points that are substantially farther away from 
the linear relation than expected. For a non-thermal population 
distribution, a different rotational temperature is anticipated for 
the J - 13i-13()A" line and thus the LTE model is not applicable 
to all the data. In addition, the observations were not simultane- 
ous and variability in the outgassing rate is likely to be present. 
If the J - 13i-13oA line is excluded from the rotation diagram 
analysis, a rotational temperature of 40 ± 3 K is derived (dashed 
line in Fig. [3]). 

Rotational temperatures of 76 + 2 K in HCN and 93 + 2 
K in HjO were obtained on 19 January 2005 with the Near 
InfraRed echelle SPECtrograph (NIRS PEC) instrument at the 
Keck II telescope (Bonev et al. 2009| . Kawakita & Kobayashi 
( 2009| l found a rotational temperature of 85 + 5 K in HjO and 



90+7 K in CH4 on 30 January 2005. Retrievals of the ortho-to 
para ratios in H,0 in C/2004 Q2 (Machholz) at infrared wave- 
lengths are consistent at 95% confidenc e level with r^pin > 34 K 



(IBo riev et aL]|2007 ) and T.pin > 27 K ( [Kawakita & Kobayashi" 
2009|l. From the A/E/F relative abundances in CH4 the lower 



limits are Tspin > 35 - 3 8 K (Bon ev etal.|2009[ ) and T^pin > 27 K 



( [Kawakita & Kobayashi 2009). I'he nuclear spin temperature 
is defined as the rotational temperature retrieved from a given 
spin isomers abundance ratio assuming their respective popu- 
lation distributions are in LTE conditions. These observations 
at infrared wavelengths prove the collision-dominated inner re- 
gion of the coma where the population levels are roughly in ther- 
mal equilibrium, and the rotational temperature is expected to be 
similar to the kinetic temperature only in the collisional region 
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( |Bockelee-Morvan et al.|1994a[ ). However, the derived rotational 
temperature of CH3OH at 338 GHz is likely to be different from 
the kinetic temperature because of relaxation of the population 
levels toward fluorescence equilibrium. We could not determine 
the rotational temperature for other gas species because no other 
simultaneous line observations were obtained. Therefore a ki- 
netic temperature of 60 K was assumed in the excitation model- 
ing for other species. 

3.2.2. Production rate from the rotational diagram 

We derived a CH3OH beam-averaged column density of 
{N}/Z{T,od = (2.22 + 0.45) x 10" cm^^ from the intercept 
of the linear fit of the J - 7-6 transitions in the rotational di- 
agram at = 0. Assuming that the population of methanol 
levels is close to a Boltzmann distribution, the derived col- 
umn density is a good approximation of the real distribution 
in the coma. Accordingly, we omitted the transition J - 13 1- 
13()A' in deriving the column density for the reasons outlined in 
Sect. 



3.2.1 



We computed the partition function evaluated at the 
rotation temperature by performing a logarithmic interpolation 
between the values from |Villanueva et al.| ( |2012) l. These parti- 
tion functions were computed using two million CH3OH levels 
for the range T = 1-500 K. The resulting column density is 
(A^) - (8.8 ± 3.2) X lO'-' cm"^ including the 1-cr uncertainty 
in the partition function introduced by Ty^t- A methanol produc- 
tion rate of 2ch,oh = (6.8 + 2.5) X 10^^ molec. s"' was derived 
from the column density assuming a Haser parent molecule dis- 
tribution with a constant expansion velocity of 0.75 km s"' and a 
beam-pointing offset of 8" from the nucleus, corresponding to a 
mixing ratio relative to water of ~ 2.5%. On the other hand, from 
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Fig. 3. Rotation diagram for CH3OH lines in comet C/2004 Q2 
(Machholz) including 1-cr uncertainties. The column density of 
the upper level divided by its degeneracy in logarithmic scale is 
plotted against the energy of the upper level. Filled (blue) cir- 
cles are A^ and A -CH3OH lines, empty (red) circles denote E- 
CH3OH transitions and the empty (red) triangle indicates the 3-cr 
upper limit on the J = 74-64E transition at 338.530 GHz. The 
solid line shows the best linear fit with a derived rotational tem- 
perature of 85 + 7 K. The dashed line shows the best linear fit 
excluding the J = 13i-13oA" transition at 342.730 GHz (energy 
of the upper state of ~ 228 K relative to the ground state) with a 
derived rotational temperature of 40 ± 3 K. 
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Fig. 4. HCN (4-3) multiplet at 354.505 GHz observed with the 
CTS on 12.96 UT January with a 6420 s total on-source integra- 
tion time. The vertical axis is the main beam brightness temper- 
ature scaled by the beam efficiency of SMT and the horizontal 
axis is the Doppler velocity in the rest frame of the nucleus. 



3.1 a production rate of 



the excitation model described in Sect. 

GcHiOH = (5.5 + 0.6) X 10^^ molec. s"' was estimated from the / 
- 7-6 transitions assuming that Tkin - Trot and the same point- 
ing offset, which corresponds to a relative abundance relative to 
water of ~ 2.1%. This mixing ratio represents an intermediate 
value in comparison with other Oort Cloud comets. 



3.3. HCN 

The HCN transitions are the brightest emission lines in 
cometary atmospheres for ground-based submillimeter obser- 
vations fSc hloerb et al.||1987[ |Bockelee-Morvan et all|1994b[ 
|Drahus et al.||2010 ). Its production rate has been found to be 
in a roughly constant ratio with respect to water, therefore it 
is often used to derive relative molecular abundances. Figure |4] 
shows the HCN J - 4-3 rotational line at 354.505 GHz ob- 
served by the CTS. The line is detected with a very high S/N 
of ~ 100 after about two hours of integration. The intensity and 
velocity shift of the observed HCN emission line intensity are 
shown in Table [3] From the width of the line we obtain an esti- 
mate of 0.75 km s"' for the coma expansion velocity. Expansion 
velocities in the range Vexp - 0.5 - 0.8 kms"' are typical for 
comets at this heliocentric distance and for this gas production 



rate (Tseng et al. 2007). Because HCN has a short rotational 



lifetime, it is important to model its excitation including col- 
lisional and radiative processes, which control the population 
levels in the outer coma. We derive an HCN production rate of 
(2.26 ± 0.02) X 10^^ molec. s"' assuming the kinetic tempera- 
ture to be equal to 60 K. The water production rate was assumed 
to be 2h,o ~ 2.5 x 10^^ molec. s"', which corresponds to an 
HCN/HjO mixing ratio of about 0.08%. This value is moder- 
ately lower than the standard ratio of 0.1% observed at radio 
wavelengths for several comets over a wide range of heliocen- 
tric distances ( |Biver et al.|2002t[Bockelee-Morvan et al.|2004[ i. 
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Fig. 5. Averaged spectrum of the H'-'CN (4-3) multiplet at 
345.340 GHz observed on 13.20 and 16.16 UT January with 
the low-resolution AOS. The vertical axis is the calibrated main 
beam brightness temperature and the horizontal axis is the 
Doppler velocity in the comet rest frame. This line is outside 
the CTS and AOSC frequency ranges. 



3.4. H^^CN 

Owing to the strong activity of C/2004 Q2 (Machholz) and 
its close approach to Earth during our observational campaign, 
it was possible to search for weak transitions of rare isotopic 
species. The H'^CN J = 4-3 transition at 345.340 GHz was 
observed simultaneously with the CS and CO transitions on 
14.20 and 16.16 UT January with a total integration time of 
~ 640 minutes, and is marginally detected by averaging the 
observations from both days (see Fig. [5]|. The inferred pro- 
duction rate is (3 ± 1) x 10^"* molec. s"'. The derived ^^C/"C 
isotopic ratio of 97 + 30 is consistent with the standard so- 
lar value of 89 observed in various comets. Isotopic fractiona- 
tion is very sensitive to chemical and physical conditions where 
the molecules condensed and therefore provide important clues 
about how cometary material formed in the early solar system. 
If comets originate from different regions in the solar nebula, 
they could display variation in their isotopic composition de- 
pending on the local temperature. Additional measurements of 
the '^C/'^C ratio i n HCN have been obtained in comets C/ 1995 
Ol (Hale-Bopp) Pewitt et al.||1997| JZiurys et al.|[T999]l and 
17P/Holmes ( Boc kelee-Morvan et al.|[2008| l, which are consis- 
tent with a solar value. The carbon isotopic composition of sev- 
eral comets of different dynamical families also agrees with a 



solar value from measurements of C, and CH (see Jehin et al. 



2009 Bockelee-Morvan 201 1, for reviews on the isotopic com- 
position of cometary volatiles). 

3.5. HNC 

Hydrogen isocyanide (HNC) is a metastable isomer of HCN and 
has been detected in several comets with a measured HNC/HCN 
ratio ranging between 0.03-0.3 ( Lis et al.|2 008 ). Figure|6]shows 
the HNC spectrum in C/2004 Q2 (Machholz) observed by the 
CTS. A production rate of (7 + 2) x 10^"* molec. s"' is derived, 
corresponding to a relative abundance of 3.1% relative to HCN, 
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Fig. 6. Averaged spectrum of the HNC (4-3) multiplet at 
362.630 GHz observed on 14.17 and 15.14 UT January with the 
CTS. The vertical axis is the main beam brightness temperature 
and the horizontal axis is the Doppler velocity in the comet rest 
frame. The effective resolution after smoothing is 200 m s"' . 



a typical abundance for moderately active comets. The origin 
of cometary HNC is still a debated topic. HNC is an abundant 
species that may be produced by chemical reactions from par- 
ent volatile species in the coma. Models have shown that HNC 
could also be formed in the coma via ion-neutral and isomeriza- 
tion chemical reactions (Irvine et al. 1998 ). Consequently, deter- 
mining the production rate is important to distinguish between 
single parent volatile and multiple precursors for HNC. Since 
we do not have precise estimates on its distribution, we assumed 
the Haser formula for a parent molecule in our computation of 
the production rate. 



3.6. 



H2CO 



Formaldehyde is a common molecule in the interstellar medium 
and has been detected in more than 20 comets at radio wave- 
lengths since the first unequivocal detection at millimeter wave- 
lengths o f the 3i2— 2ii line at 226 GHz in comet C/1989 XI 
(Austin) ( |Colom et al. 1992| l. Solid-phase hydrogenation reac- 
tions on the surface of the nucleus could account for the pres- 
ence of formaldehyde and have been confirmed recently by lab- 
oratory experiments ( Nagaoka et al.| |2005; Hidaka et al. 2009)l. 
Figure [7] shows the HjCO 5 15^14 transition at 351.769 GHz in 
comet C/2004 Q2 (Machholz) obtained by the AOS. CTS and 
high-resolution AOS data are not available for this transition. It 
is possible that the frequency scale in this range of the AOS can 
be slightly inaccurate. For the acousto-optical backends, the cor- 
respondence between channel number and frequency is not quite 
linear, and the correction may be important and ill-known on the 
border of the band. This line is detected in the lower sideband of 
the AOS at about 300 MHz from the central frequency with a to- 
tal bandwidth of ~ 1 GHz. The same HjCO line was observed at 
the Caltech Submillimeter Observatory (CSO) on similar dates 
(12.1 and 13.1 January) where the velocity shift toward the blue 
wing was found to be about twice smaller (Biver et al. in prepa- 
ration). 

HjCO has ortho- and para spin isomers that differ in the 
alignment of the nuclear spin of the H atoms (parallel in the or- 
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Fig. 7. HjCO (5 1 5^ 14) transition at 351.769 GHz observed on 
12.96 UT January with the lower resolution AOS. The vertical 
axis is the main beam brightness temperature and the horizontal 
axis is the Doppler velocity in the comet rest frame. This line 
is outside the CTS and AOSC frequency ranges. The frequency 
scale in the range close to the border of the AOS can be slightly 
inaccurate. 



tho state and anti-parallel in the para state). We assumed that 
cometary formaldehyde has an ortho-to-para ratio of 3 accord- 
ing to the nuclear spin statistical weights ratio. HjCO is a short- 
lived molecule and its spatial distribution has been found to dif- 
fer from that expected for a parent molecule. An extended source 



of H^CO has been identified in several comets (Meier et al. 1993 



Bockelee-Morvan et al.|1994m. Assuming that HjCO originates 



either from the nucleus or from an extended source with a par- 
ent scale length of Lp - 8000 km, production rates vary between 
(3.87+0.20)xl026molec. s"' and (1.09+0.06)xl02^ molec. s"', 
respectively, as inferred by a model with an 8" average point- 
ing offset (Biver 1997 ). The first value represents a lower limit 
because it is expected that direct release from the nucleus un- 
derestimates the real H2CO production rate. The corresponding 
Gh^co/GhjO mixing ratios are 0.14% and 0.41% for the parent 
source and extended HjCO source production. These values are 
intermediate in comparison with the observed mixing ratios in 
other Oort Cloud comets ( |Bockelee-Morvan|201 l| l. 

3.7. CO 

We show the carbon monoxide J - 3-2 line in Fig. [8] The line 
profile is strongly asymmetric with a peak in the redshifted wing. 
There is some observational evidence that the CO lines may be 
affected by self-absorption effects in CO-rich and very active 



comets, which have been included in our analysis (Bockelee- 
|Morvan et al.||20l"0l l. For the observed CO transition in comet 
C/2004 Q2 (Machholz) our model predicts an optical depth of 
T ~ 0.2. Therefore, the observed excess emission at redshifted 
frequencies is attributable to kinematic effects by the presence of 
a jet roughly in the antisolar direction since the solar phase an- 
gle was small during the observations. An anisotropic outgassing 
radiative transfer model with a rotating jet-like structure would 
provide a better fit to the line shape (see e.g. Bockelee-Morvan 
et al.|200^ I. Based on an isotropic coma distribution and using 
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Fig. 8. CO (3-2) transition at 345.796 GHz observed on 14.94 
UT January with the CTS. The asymmetric line profile suggests 
a strong anisotropic outgassing with preferential production in 
the direction opposite to the Sun. 



rate of (7.0 + 0.6) x lO"^' molec. s"' . CO is thought to be a parent 
molecule with a small additional contribution from dissociation 
of CO2 and we used the Haser parent molecule distribution to 
derive the CO production for simplicity. 

Oort Cloud comets are expected to be enriched in hyper- 
volatile species such as CO and the abundance ratio of HjCO 
and CO has been observed to be correlated in several comets. We 
find a relative abundance of CO with respect to HCN of about 
31 and a mixing ratio of ~ 2.6% relative to HjO in C/2004 Q2 
(Machholz) . Observations of the hypervolatile CO in comets 
at infrared and radio frequencies show a wide range of abun- 
dances ( |Bockelee-Morvan et al.|2004l l. CO abundances observed 
at radio wavelengths range between ~ 0.4-20% relative to HjO 



a kinetic gas temperature of 60 K, we deduce a CO production 



(Bockelee-Morvan |201 1 ). Generally, these lines are relatively 
weak because of the small CO dipole, although they are the 
most readily observed cometary lines at large heliocentric dis- 
tances. The CO abundance in eight Oort Cloud comets observed 
by infrared ground-based spectroscopy ranged between ~ 1- 
20% relative to water (Mumma et al. 2003 ). Thus the measured 
2co/2h,o mixing ratio lies toward the low end of the observed 
range in Oort Cloud comets. Depletion of CO would therefore 
suggest that C/2004 Q2 (Machholz) has been injected into the 
inner solar system before and is not a dynamically new comet. 

3.8. CS 

The CS radical is believed to be produced in cometary coma 
from carbon disulfide (CSj) and has been observed in several 
comets in the UV and radio wavelengths ( |Feldman et al.|2004t 
Bockelee-Morvan et al.|20()4 i. Figure |9] shows the observed CS 
J - 7-6 emission line at 342.883 GHz. The line profile is asym- 
metric with a peak in the redshifted side that resembles the CO 
J = 3-2 line shape shown in Fig.|8] This suggests an anisotropic 
outgassing in the anti-sunward direction. It is expected that the 
line is optically thin and absorption effects in the foreground are 
negligible according to our model. To obtain the density profile 
we used a photodissociation rate of yScs = 2.5 x lO^^molec. s~' 
at Th = 1 AU derived from spectroscopic observations ( |Boissier| 
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Fig. 9. CS (7-6) emission line at 342.883 GHz observed on 
14.94 UT January with the CTS. A double-peaked profile is ob- 
served. The velocity scale in the horizontal axis is in the comet 
rest frame. 



et al.|2007 [Biver et al.|2011) l. This photodissociation rate corre- 
sponds to a scale-length of 28 000 km, assuming an expansion 
velocity of 0.75 kms ' obtained from the width of the HCN 
line. A CS production rate of (1.15 + 0.04) x lO^*" molec. s ' 
is estimated from a spherically symmetric model with direct re- 
lease from the nucleus. Assuming a distributed source in the 
coma with scale lengt h - 650 km, w hich corresponds to 
production from CSj ( Boissier et al.|20()7 1, a production rate of 



(1.29 ± 0.04) X 10^^ molec. s'' is inferred. The derived mixing 
ratio with respect to water of ~ 0.04% is lower than the typical 
value measured in other comets at radio wavelengths ( |Bockel6e^ 
|Morvan|201 1] [Crovisier et al.|2009a| l. 



3.9. Short-term outgassing variability 



Several comets have shown periodic variations in the outgassing 
that are associated with the rotation of the nucleus, combined 
with a long-term seasonal variation with a peak emission close 
to perihelion. In some cases, the rotation periods of cometary 
nuclei have been determined by measuring this periodic out- 
gassing variability using various techniques (see [Samarasinha 
et al. 2004| for a review). Rotation periods have been derived 
from the periodic variability of the HCN production rate in 
comets 9P/Te mpel 1 (Biver et al. ,2007a) , 73P-C/Schwassmann- 
Wach mann 3 ( |Drahus et al.||201()| l and 2P/Encke Pockers et al.| 
201 l| l. To address the question of short-term periodic variability 
in the production rates of C/2004 Q2 (Machholz), we selected 
the HCN and CS observations from one night, 12-13 January, 
which have sufficient S/N. 

The production rates for each individual scan of the HCN 
observations are shown in Fig. 10 The HCN J = 4-3 line 



was detected in each single 8-minute scan with an S/N of ~ 
15. We observe an increase of the production rate by roughly 
40% during the 2-hour observing interval with a mean value 
of (2.09 + 0.02) X 10^'' molec. s"', although the first scans at 
the beginning of the observations were obtained at lower eleva- 
tion starting at 45° and could be affected by stability problems 
in the instrument or changes in the main beam efficiency. The 
HCN evolution suggests a periodic variation in the activity of 
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Fig. 10. HCN production rates in comet C/2004 Q2 (Machholz) 
as a function of time on 12-13 January 2005. Each data point 
represents an 8-minute observation. 



the comet potentially induced by the rotation of the nucleus due 
to the presence of active regions on the surface. The large am- 
plitude of the variation can be explained by an intrinsic activity 
change as one or more active areas are exposed to solar radia- 
tion during the rotation cycle. This variation is compatible with 
the rotation periods of 9.1 ± 0.2 hours derived by Reyniers et al. 
( |2009| l from broadband filter optical photometry, 9.1 + 1.9 hours 
reported by Sastri et ar] ('2O05) from studying dust fans visible in 
/?-band images, 17.60 + 0.05 hours from narrowband CN images 
reported by Farnham et al.| ([2007 ), and 17.8 + 0.5 hours by deter- 
mining the itmer coma morphology from photometric variations 
( [Manz ini et al. |2012| l. However, we are not able to perform a de- 
tailed perioJicity analysis of the HCN production curve due to 
the insufficient temporal coverage because the total time span of 
our observations is less than one night. Hence it is not possible 
to determine if the brightness increase is aperiodic phenomenon 
or due to non-periodic activity changes. 

Figure 1 1 shows the evolution of the CS outgassing rate on 
13 January during a period of 7.5 hours. The CS line is rela- 
tively strong and has one of the highest peak antenna temper- 
atures after the HCN line. CS is detected in most of the indi- 
vidual 8-minute scans with an S/N between 3-6. The comet dis- 
played a strong variability with a 40% variation around the mean 
value. It is likely that this brightness variability is not a periodic 
event. Small outbursts or activity changes are expected as comets 
approach the Sun, which could explain the observed variability 
(A'Hearn et al.^2005). It is also possible that the strong fluctua- 



tions in the CS production could be partly explained by pointing 
uncertainties that affect the derived production rates combined 
with instrument instabilities. Observations of the CS line on 16 
January have a lower S/N and only in a fraction of the individual 
scans the Une is detected with at least a 3-cr confidence level. 



4. Discussion 

We have observed several molecular species in comet C/2004 Q2 
(Machholz) when it was close to perihelion in the sub-millimeter 
wavelength range with the SMT during six consecutive nights 
in January 2005. These observations led to the detection of sev- 
eral CH3OH, HCN, H'^CN, HNC, HjCO, CO, and CS rotational 
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Fig. 11. Time evolution of the retrieved CS production rates in 
comet C/2004 Q2 (Machholz) on 13 January 2005. Each data 
point represents an 8-minute observation with statistical uncer- 
tainties. 



lines. Our main goal was to measure the relative production rates 
of several parent and daughter volatiles in this comet and esti- 
mate the CH3OH rotational temperature. We derive an HCN pro- 
duction rate of (2.26 ±0.03) x 10^^ molec. s"' at heliocentric dis- 
tance of 1.2 AU, which corresponds to a 2hcn/6h,o ~ 0.08%, 
using a spherically symmetric radiative transfer numerical code 
that includes collisional effects between neutrals and electrons 
and radiative pumping of the fundamental vibrational levels by 
solar radiation ( |Bockelee-Morvan & Crovisier||1989| ). In addi- 
tion to the statistical noise, the precision of the production rates 
is affected by the limited pointing accuracy at the time the obser- 
vations were performed. An 8" pointing offset was included in 
the computation of the production rates, which results in ~ 20% 
uncertainty in the derived values. Mixing ratios relative to hy- 
drogen cyanide and water are listed in Table|2]for the production 
rates derived from the weighted average value of the observa- 
tions on different dates during our campaign. The asymmetric 
shape of the CO and CS line profiles suggests that there is pref- 
erential outgassing from the anti-sunward side of the nucleus, 
while the HCN and HjCO line profiles are fairly symmetric, al- 
though the latter is shifted toward the blue wing perhaps due to 
an instrumental effect. Using the rotational diagram technique, 
we retrieve a cold rotational temperature of 40 ± 3 K for the 
methanol energy levels sampled by the J = 7-6 transitions. 
Retrieved molecular abundances relative to water are com- 



parable to those obtained in Oort Cloud comets ( Biver et al. 
[2OO2I ICrovisier et"ar]|2009b| |Dello Russo et al.|]20TTp . Our 
Ghcn/2h,o mixing ratio of ~ 0.084% is slightly lower than 
the typical value of 0.1% observed at radio wavelengths, and 
50% lower than those found by |Bonev et al.| ( |2009| l from av- 
eraged 28 November 2004 and 19 January 2005 observations 
and by (Kobayashi & Kawakita 2009) from observations on 30 
January 2005 using the same instrument. The observed differ- 
ence between our measurement of the HCN mixing ratio and 
those derived at infrared wavelengths is fairly typical. The mix- 
ing ratio of CH3OH relative to HjO agrees with the infrared 
measurement in comet C/2004 Q2 (Machholz) by Bonev et al. 
p009| l performed on 19 January 2005 at rh = 1.208 AU that is 
closer in time with our observations. Considering the pointing 



uncertainty that introduces an error of ~ 20%, our measurement 
is consistent within 1-cr with the revised abundances retrieved 
from the infrared observations on 19 and 30 January 2005 by 
jVillanueva et al. (2012) using a new line-by-line model for the 
V3 fundamental band of CH3OH. The averaged CO production 
rate measured on 13 and 16 January agrees with that obtained 
by IBonev et"aLl p009| l on 29 November 2004 at rh = 1.493 
AU ((6.3 + 0.3) X 10^ molec. s ') within confidence limits. In 
contrast, the Qco/Qh^o mixing ratio of ~ 2.6% is almost a 



factor of two lower than that derived by Bonev et al. (2009 1. 
For a parent molecule distribution, the derived Qh^co/Qh^o ra- 
tio is 0.14%, and 2h,co/6hcn is 1.7%, which is more reliable 
since these molecules were observed simultaneously. These val- 
ues are intermediate between those measured by Bonev et al. 
( |2009| and [Kobayashi & Kawi^ita| ( |2009| . On the other hand, 
the values inferred from a daughter molecule extended source 
distribution with a scale length of Lp = 8000 km are about a 
factor of 3 higher. This illustrates that the derived production 
rates depend strongly on the assumed scale length of the parent 
molecule. This is particularly the case for HjCO where the beam 
size is slightly smaller than Lp, so it may be possible that this 
value overestimates the production rate. An HNC abundance of 
3.1% with respect to HCN is found assuming direct release from 
the nucleus with a Haser distribution (3.0% if infrared pumping 
of the fundamental vibrational levels is not considered), which 
suggests that this molecule may be destroyed by chemical re- 
actions. This value is compatible with those observed in other 
comets given the dependence on heliocentric distance - except 
73P/Schwassmann-Wachmann and the very active Hale-Bopp 
(C/1995 Ol) (see Lis et al.|2 008). 

Outgassing variations induced by the nucleus rotation are ex- 
pected to appear from non-sphericity of the nucleus or the pres- 
ence of active region areas on the surface. The variability ob- 
served in the HCN and CS outgassing rates on 12-13 January 
is affected by inaccurate pointing, which introduces an uncer- 
tainty of about the expected variation caused by the rotation of 
the nucleus. HCN production rates show a uniform brightness in- 
crease over a period of two hours, which could be a periodic phe- 
nomenon consistent with the rotation period of the nucleus de- 
rived using different methods (S astri et al.[ 2005, Farnha m et al.| 
5007t [Reyniers et"aLp009, .Manzini et al.|2012) , but our obser- 
vations do not provide sufficient phase coverage to constrain the 
rotation period. 

Comets are the most pristine objects in the solar system and 
have not undergone substantial thermal processing. Cometary 
ices are more sensitive to thermal processing than dust. Hence 
their abundances provide indications about the formation and 
evolution of material in the early solar nebula. Dynamically new 
comets are expected to be enriched in volatile species, and CO is 
the most volatile component observed in C/2004 Q2 (Machholz). 
However, our observations show a relative depletion of CO and 
and an intermediate-range mixing ratio of HjCO and CH^OH 
compared to the range of values measured in other comets ( Biver 
et al.||20(j2 l. It is believed that the composition of icy material 



on the surface of the nucleus is altered by the exposure to so- 
lar radiation. Thus, these observations provide a hint about its 
thermal/dynamical history and suggest that the comet has visited 
the inner solar system previously. Formation regions of comets 
of Oort Cloud comets vary from 5-30 AU from the protosun 



according to the standard "Nice model" scenario (Morbidelli 
|et al. '200¥) and can be constrained by observations of the chem- 
ical composition. The relative depletion of several volatiles in 
C/2004 Q2 (Machholz) agrees with the determination of the for- 
mation region in the inner region of the solar nebula compared 
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with other Oort Cloud comets, derived from the lower limit of 
the retrieved nuclear spin temperatures and dynamical models 
of the evolution of planetesimals in the solar system ( ^Kobayashi] 
[&Kawakita 2009; Kawakit a & Kobayashi|20"09l ). 

Submillimeter spectroscopy of cometary atmospheres is a 
useful approach for studying the diversity of species that sub- 
limate when a comet approaches the Sun. These observations 
demonstrate the capabilities of the new CTS based on digital 
technology installed at the SMT to support solar system obser- 
vation programs, further established with the cometary observa- 



tions presented in Kiippers et al. ( 2004 1; Villanueva & Hartogh 
^2O0^■ [Drahus et al.| ( |2010| l; |Pagamni et al.| ( |2010| l; [Jockers 



et al.|(2011|l, and observations of the Venusian mesosphere in 
Rengel et at] ( |2008a|b[ ). Our observations of comet C/2004 Q2 



(Machholz) were analyzed for noise properties, standing wave 
and spectral line properties as compared to model predictions. 
The spectrometer is found to perform according to its original 
specification in terms of sensitivity, spectral resolution and sta- 
bility. 
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Table 3. Rotational emission lines in comet C/2004 Q2 (Machholz) observed by SMT on 13-16 January 2005. 



Molecule Label" Date* 


Transition Frequency'' 


J T^B dv 


Velocity shift'' 
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(UT) 
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(Kkms-') 


(m s"') 
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HCN" 
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4-3 


354.505476 


3.962 + 


0.033 


-31 + 


7 


(1.87 


+ 0.02) X 10'" 


(2.26 + 


0.02) X 10'" 


0.34 X 10'" 


H2CO' 
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5l5^14 


351.768645 
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0.035 
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(9.64 


+ 0.50) X 10-*^ 
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Notes. Observed flux densities integrated over velocity, Doppler velocity shifts and production rates derived using the comet model described in Sect. |3.1 | are shown with statistical uncertainties. 

Label s in Fig. [T]**' Mid- time of the observations recorded as fractional days in UT. ''' The line frequencies were obtained from the latest online edition of the JPL Molecular Spectroscopy 
Catalog j Pickett eTal.|l998" . The velocity offsets are computed with respect to the optocenters of the complete components for blended lines. Production rates assuming a pointing offset of 2". 

Production rates assuming a pointing offset of 8". The uncertainty in the production rates are 1-cr uncertainties considering the receiver sideband gain ratio, main beam efficiency and kinetic 
temperature errors. Line intensities are the sum of the hyperfine components. Production rates were derived assuming a daughter product extended source distribution with a scale length of 
Lp = 8000 km for H2CO and Lp = 650 km for CS. Weighted average of observations obtained on 13.20 and 16.16 UT January. 



